Generating and Adding Flows on Locally Complete 

Metric Spaces* 



Hwa Kil Kim and Nader Masmoudi 
March 8, 2013 

Abstract 

As a generalization of a vector field on a manifold, the notion of an arc field on 
a locally complete metric space was introduced in [4]. In that paper, the authors 
proved an analogue of the Cauchy-Lipschitz Theorem i.e they showed the existence 
and uniqueness of solution curves for a time independent arc field. In this paper, 
we extend the result to the time dependent case, namely we show the existence and 
uniqueness of solution curves for a time dependent arc field. We also introduce the 
notion of the sum of two time dependent arc fields and show existence and uniqueness 
of solution curves for this sum. 



1 Introduction 

Vector fields play an important role on manifolds. In particular they allow the study of 
dynamics on the manifold. On metric spaces and in the absence of a differential structure, 
the notion of arc fields was introduced in [3]. Under some regularity assumption, the authors 
of [3] proved the existence of solution curves for a time independent arc field. Their result 
can be seen as an extension of the Cauchy-Lipschitz Theorem. The goal of this paper is to 
define a notion of sum of two arc fields and construct a unique solution curve for this sum. 
We also generalize [1] to the time dependent case. 

Let us also mention that the generalization of the notion of differential equations from 
manifolds to metric spaces is a natural question. In this direction, there are many other 
approaches which can be found in [3], [10], [12], [13], [13] and [5]. A basic idea that all 
approaches have in common is to replace the concept of a vector field by a suitable family of 
curves (herein called an 'arc field' following [1]) each of which supplies the direction of travel 
at the point from which it issues. We borrow the idea of [3] which shows the existence of 
flows corresponding time independent arc fields on locally complete metric spaces whereas 
all others have predominantly assumed that the underlying metric space is locally compact. 
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Let us now explain our motivation behind this work. In [9], we study systems coupling 
fluids and polymers. In it most generality the phase space for the polymers is given by a 
metric space (see [8]). When the phase space of the polymers is a manifold, we get a system 
coupling the Navier-Stokes equation for the fluid velocity with a Fokker-Planck equation 
describing the evolution of the polymer density (see for instance [HI [7]). The coupling comes 
from an extra stress term in the fluid equation due to the polymers. There is also a drift 
term in the Fokker-Planck equation that depends on the spatial gradient of fluid velocity. It 
can be seen that the Fokker-Planck equation has a flow structure on the set of probability 
densities of polymers. More specifically, let Ai be the set of all Borel probability measures 
defined on the manifold (phase space of polymers) then we can put a metric structure on 
Ai using the Wasserstein distance. Once Ai is equipped with the Wasserstein distance, the 
Fokker-Planck equation can be considered as the sum of two flows on Ai. One is the gradient 
flow corresponding to the entropy functional on Ai and the other one is a drift term which 
is generated by the spatial gradient of the fluid velocity which depends on time. If the phase 
space of polymers is not a manifold but just a metric space then we don't have Fokker-Planck 
equation any more. But the flow interpretation is still available to describe the evolution of 
the polymer density if we know how to generate and add flows on metric spaces. Achieving 
this is one of the goal of this paper. 

We briefly summarize the contents of each section. In section 2, we study time dependent 
arc fields, solution curves, and sufficient conditions under which we can prove the existence 
of solution curves for arc fields. We also show the continuous dependence of solutions on 
initial conditions from which we can get the uniqueness of the solution curve. In section 3, 
we introduce the notion of solution curve for the sum of two arc fields. By imposing a kind 
of commutation law on two time dependent arc fields, we prove the existence of solution 
curves. We also get the uniqueness of a solution curve to the sum of two arc field by showing 
the continuous dependence of solution curves on the initial conditions. 

2 Generating Flows 

2.1 Time dependent arc fields 

Let X be a locally complete metric space with a metric d. 

Definition 2.1. A time dependent arc field on X is a family of maps $(■; •, •) : [0, oo) x X x 
[0, 1] — > X such that for all t G [0, oo), x G X, we have $(£; x, 0) = x, 

d{${t;x,h),®{t;x,k)) 

p(x, t) : = sup — < oo 

h+k \h-k\ 

and the function p(x, t) is locally bounded, namely for all t, x G [0, oo) x X, there exist 
r, / > such that 

p(x,t;r,l) := sup {p(y, s)} < oo. 

yeB(x,r),\t-s\<l 
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One can interpret <3>(i; b, •) : [0, 1] — > X as a curve on A starting from b to $(£; 6, 1). This 
gives the direction of the curve in some sense. Notice that, for fixed b G X, the direction 
given by <3>(i; b, .) depends on the time i. Besides, p(x,t) can be understood as the upper 
bound on the speed of the curve •). For the convenience, we will use the notation 

:=<&(*; 6, h). 

Definition 2.2. For given a £ X and £ G [0, oo), a solution curve of $ with initial position 
a at time t is a map a : [t, £ + c) — >■ X (for some c > 0) such that a(t) = a and for each 
s G [£, t + c) 

lim *(£±^1IM£))) = o p.!) 
h^o+ h 

We introduce some conditions on the time dependent arc field $. Motivations for Con- 
dition A and B were already given in [I]. Condition C is about the time regularity of $. 

Condition A: There is a function A : X x X x [0, 1] — > (—00,00) such that for each 
a G X and t G [0, 00), there are constants r a > 0, e a G (0,1] and T t > t such that A is 
bounded above on B(a,r a ) x B(a,r a ) x [0, e Q ] and 

d($£(ai),$£(a 2 )) <d(a 1 ,a 2 )(l + /iA(a 1) a 2 ,/ i )) (2.2) 

for all ai,a 2 G B(a,r a ), h G [0,e a ] and s G [t,T t ]. 

Condition B: There is a function Q : X x [0, 1] x [0, 1] — >■ [0, 00) such that for each 
a G X and t G [0, 00), there are constants r a > 0, e a G (0, 1] and T t > t for which f2 is 
bounded on B(a,r a ) x [0, e a ] x [0, e Q ] and 

d(<f> s l+h (b), n o $f (6)) < hg(h hMb, /, ^) (2.3) 

for all b G B(a, r a ) , l,h G [0, e ] and s G [t, T t ] where (7 : [0, e ] x [0, e„] — )■ [0, 00) satisfies 

lim g(l,h) = Q and p(2 _i ,2 _i ) < 00 (2.4) 

Condition C: For each a G A and t G [0, 00), there are constants r a > 0, e a G (0, 1], 
T t > t, < a < 1 and C > 0, such that 

d(^(b),^(b))<Ch\ Sl -s 2 \ a (2.5) 

for all b G B(a, r a ) , h £ [0, e ] and s±, s 2 G [t, T t ]. 

Remark 2.3. Once we have fixed a G A, t G [0, 00) and fixed constants r a , e a , T t then functions 
A and fl are bounded above. We denote upper bound of A (respectively Q) by Ka(Kb)- 

As a simple observation, by combining Condition B and C, if b, G B(a,r a ) then we 
have 

rf($? + ,(6), & h +l o $•(&)) < d(^(6), n o <&? (6)) + d(& h o $f(6), o $•(&)) 

< hg{l,h)K B + Chl a 

= h(g(l,h)K B + Cl a )=:h~g(l 1 h) (2.6) 
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Lemma 2.4. For a given a G X andt G [0, oo), letr a , e a andT t be the constants in Condition 
A, B and C. Ifb u b 2 G B(a,r a ), 2h G [0,ej, s, s + h G [t,T t ] and & h (bi), $ s h (b 2 ) G B(a,r a ) 
then we have 



d{% +h o $*( fel ), & 2h (b 2 )) < d(b x , 6 2 )(1 + WT A ) a + ^(a, fc) 



(2.7) 



b*(A) 



Figure 1: 

Proof. Triangle inequality gives 

d(K +h o *'M, < d(<S>* h +h o n(bi),n +h o n(b 2 )) + d($ s 2h (b 2 ), <S>° h +h o $*(& 2 )) 

(2.1 



For the first term in the right hand side of ( 12. 8ft . we use Condition A twice 
d« h o 9' h (h), n +h ° < d(*h(h), nihMl + WO) < d(h, b 2 )(l + Wf x ) a (2.9) 



For the second term, we exploit (12.61) to get 

d($ s 2h (b 2 ),®l +h o<t>l(b 2 ))<h~g(s,h) 



(2.10) 



We combine (12TB]) . (T2~9|) and (12~T0D to finish the proof. 



□ 
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Remark 2.5. In general, we have 

d(^ o ^(60, <^ 1+ , 2 (6 2 )) < d(h, 62) (1 + ^^)(1 + Zi 2 A A ) + /^(/n, fc 2 ) (2.11) 

Lemma 2.6. For a given a G X andt G [0, 00), letr a , e a andT t be the constants in Condition 
A, B and C. Ifbi,b 2 G B(a,r a ), h, I + h G [0, e a ], s, s + Z G [t,T t ] and $f(o 2 ) G B(a,r a ) then 
we have 

d(*f '(&i), $ s l+h (b 2 )) < d(b U $'(62)) + 62, Z, Zl) 

where rj s (bi, b 2 , Z, h) := d(bi, $>i(b 2 ))KA + g(l, h). Furthermore, we notice that r] s (bi, bi, Z, h) 
converges to as d(bi, b 2 ), I, h — > 0. 

Proof. We use Condition A and (12. 6 j) to get 

d(^(6i), $^(6 2 )) < d^Cfti), ° *f (6a)) + d(* ^ ° *,'(&*), S^fa)) 
<d(& 1 ,$?(& 2 ))(l + fciCi) + 
= d(6i, $f(6a)) + h[d(h, ^{b 2 )K A + sf(/, h)\ 

and trivially n s (oi, bi,l,h) — > as d(&i, & 2 ), Z, — ^ 0. □ 

Lemma 2.7. For a given a G A and t G [0, 00) , let r a , e a and T t be the constants in 
Condition A, B and C. Ifbi,b 2 G B(a,r a ) and h G [0, e a ], s, u G [t,T t ] then we have, 

d^Kh), $«(6 2 )) < d(h, 6 2 )(1 + Zil^O + C7i|s - u\ a 

Proof. We combine Condition A and C to get 

d^Kh), $«(6 2 )) < d^Kh), n(b 2 )) + d(<s> s h (b 2 ),$t(b 2 )) 

< d(h, b 2 )(l + hK A ) + C7i|s - u\ a 

□ 

Lemma 2.8. For a given a G A and £ G [0, 00) , let r a , e a and T t be the constants in 
Condition A, B and C. Assume b±,b 2 G B(a,r a ) and s,s + h G [i, Define a polygonal 
path p(l) : [0, h] — > X starting at b\ G A as follows; p{l) : = 3>f(&i) /or < Z < ri and 
p(Z) := ^f_^*(p(ri)) for r^ <l < r i+ \ with < r\ < ■ ■ ■ < ri < r i+ i < ■ ■ ■ < — h. Then we 
have 

d(p(l), $f (6 2 )) < d{b u 6 2 ) + h max rj s (p(ri), b 2 , n, r i+1 - n) (2.12) 

l<i<fc 

for < Z < h. Furthermore, we have 



lim max rj s {p(ri), b 2 , r h r i+1 - r { ) = (2.13) 

h->0 0<n<h 
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Proof. If < I < ri then by Lemma 12. 6[ we have 

d(p(l), $f(fe 2 )) < tZ(6x, 6a) + M&i, 6a, 0, /) 
For ri < / < r 2 , we use Lemma [2.61 twice to get 

d(p(l),^(b 2 )) = d(^ r r i( P ( ri )),^ l+l _ ri (b 2 )) 

< d(p(n), $^ (62)) + (I - ri)^(p(n), 6 2 , n, / - n) 

< d(6i, b 2 ) + nr] s (bi, b 2 , 0, n) + (i - ri)%(p(ri), b 2 , n, / - 

In general, for r, < / < r i+1 , we have 

$f (62)) < d(b!,b 2 ) + rir) s (bi, b 2 , 0, n) + (r 2 - ri)r^ s (p(ri), 6 2 , n, r 2 - n) 
■■■ +{l-r i )r] s (p(r i ),b 2 ,r i ,l-r i ) 

which gives 

d(p(l), $f (62)) < d(6i, b 2 ) + h max r} s (p(ri), b 2 , r h r i+1 - n) 

l<i<k 

Equation (12.1 3p is almost trivial since d(p(ri), b 2 ),ri, r i+1 — ri converge to as h —> 0. □ 



2.2 Existence and uniqueness of a solution curve 

The proof of the next theorem is similar to the one in [2J. We can also think of it as a 
corollary of Theorem 13.61 But, to give an idea for the proof of Theorem 13.61 which is more 
complicated, we give a full proof here. 

Theorem 2.9 (Existence). Let $ : X x [0, 1] x [0, 00) be an arc field satisfying Condition A, 
B and C. For a given a G X and t G [0, 00), there exists a solution curve a : [t, t + c) — > X 
with initial position a at time t. 

Proof. For a positive integer n, we define the n-th discretized solution by 

& s (a) < s < ± 



$*+2-"(£ / J_\\ JL < s < 
^ 9_2-" V 2" / / 2" — — 



Cn(s) := < 



^ s _j. 2 -"VSnV2«// 2™ — — 2" 



i+1 



Suppose r, I > are chosen so that p(a, t; r, /) < 00. If p(a, £; r, /) = 0, then a(s) 
a solution curve. Thus we assume p(a, t; r, I) > 0, and let 



c := mm 



p(a,t;r,iy 
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a defines 



(2.14) 



It is easy to see that we have £ n (s) G B(a,r) for < s < c. This implies {^ n }^i is equi- 
Lipschitz with Lipschitz constant p(a,t;r,l). Moreover, by choosing r smaller if necessary, 
we may assume that there are constants Ka, K b and e G (0, 1] such that A(p, q, h) < Ka 
and f2(j>, h) < Kb for all p, q G B(a, r) and /', h G [0, ej. We may also assume that B(a, r) 
is a complete metric space. 




Figure 2: 

Let us first estimate the uniform distance between £ n and £ n _i. We apply Lemma 12.41 
with h = l/2 n and b\ = b<i = a to get 

2 2 1 _ 1 1 

d(£n(—),£n-i(—)) < T^fi^T^T^") 

Similarly, we apply Lemma 12.41 multiple times and get 

< 1^(1 l)h + (i + ^i) 2 l 

— 2 n 2 n 2 n 2™ 

< ^) [i + (i + §) 2 + o + §) 2 1 
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In general, for all i so that i ■ 2/2~ n < c, we have 

i-1 



3=0 

J_~(J_ 1 Jl + 2-"i^ A ) 21 - 1 

- 2 n 9 \2 n ' 2 n > (l + 2- n K A ) 2 - 1 
< l, (l + 2-^) c2 "-l 

, _ / 1 U e cJ ^ - 1 

- S'l^' 2™> K A (2 + 2~ n K A ) 



where K := £C ^^ 1 is a constant independent of n. 
So for any s G [0, c), let i be an integer such that 

2i ^ 2(i + l) 

— < s < — 

2" ~~ 2 n 

then we have 

2? 9i 9i 9? 

k-iOO) < d(Us),U ¥ )) + d(U§i),U-i(^ )) + d(^i(^),e»-i(*)) 

2z 1 1 2z 

< Li P(^) ^ - ^) + %, + Lip^-0 (S - -) 

< ±p(a,t-,r,l)+g(±,±)K (2.16) 

where we exploit the equi-Lipschitz property of £ n and (12.151) . 

Next, we exploit ( I2.16P to show that {£„} is a Cauchy sequence in the uniform topology. 
For any s G [0, c), we have 

m— 1 

d(f»(s)»£n+m(s)) < ^d(£n+j(s),£ n+i+ l(s)) 
3=0 

^ 4 , . 1 1 , 

- 2(^+3+i) P ^ a ' t;r ' ' ) +fl, 4( n +J +1 )' 2( n +J'+!)^ 

3=0 

< ^TP(a,*;^0+^ ^ 9^~ 3 , 2-^)^0 as n -)• oo (2.17) 

j=n+l 



Since £ n (<s) is in the complete space B(a 1 r), we know that £ n converges uniformly and 
we define a : [0, c) — >■ X by 

a(s) = lim £ n (s) 



It is trivial to see <r(0) = a and let us check that 

^mi±3^mm =0 (2 . 18) 

holds for all s G [0,c). Let e > and /i > be fixed such that s + h < c. From triangle 
inequality, we have 

d(v{a + h),*?>(a(s))) d(ir(s + h),Us + h)) + *Q, ffi'fofc))) 

h — h h 

, <*(*^n(*))^ + '(fr(*))) (219) 

Since £ n converges uniformly, we can choose n large enough so that 

d(g(s + h),Us + h)) d(^ a (Us)W(a(s))) e 
h h ~ 2 

We combine ( 12~T9|) and ( 12~20|) to get 



< - (2.20) 



d(a(s + h),^'(g(s))) ^ djUs + h^'jUs))) | g 
h ~ h 2 

We need to estimate the second term of (12.211) . 
Let z be such that i/2 n < s < (i + l)/2 n , then 

d(Us + h)^[ +s (Us))) < d(Us + + *0) + dftn(^ + h)M^{U^))) 

+ d{K + ^(U^)lK +S (Us))) (2.22) 

Let us estimate the righthand side of (I2.22p term by term. First, by the Lipschitz property 
of £ n , we have 

d(£ n (s + h),( n (— + h)) < (s - — )p(a, t; r a , e a ) 

< — p(a,t;r ,e a ) (2.23) 
For the second term, we use Lemma [2.81 with ^(w) := &i = &2- 

<en(^ + ^),C^(en(^)))</ io Sup^(U^ + r),en(^),r,^) (2.24) 
By using Lemma [2.71 and the Lipschitz property of £ n , we can estimate the last term 

« 2 "(^(t+ -)),^(^n(s))) < dft»(^),enW)(l + ^A) + ^ 

1 Ch 
< -p(a, i; r a , e a )(l + hK A ) + — (2.25) 
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We combine equations f l2.22p .( l2.23p . fl2.24p and (I2.25p . and assume n is large enough to have 

i ill Ch 

d(Us + h), $! +s (£n(s))) < h sup 77(e„(— + r), £ n (— ), r, — ) + — p(a, t; r a , e a )(2 + HKa) + — 

0<r<h Z Z Z Z Z 

<h sup ^ B (± + r),^),r,^) + | (2.26) 

0<r</i z Z Z Z 

We combine ( I2.2ip and (12.26|) . and let h — > with n large 

lun — ;, ; h K K JJJ < lim sup ??(£„(— + r), £„(—), r, — ) + e 

= e 

This gives fl2TT8|) . We define a : [t, £ + c) X by 

it(s) := <r(s — t), s G [t, £ + c) 

It is trivial that o~(t) = a(0) = a and (I2.18P implies that a satisfies (12. ip . This concludes the 
proof. □ 

Theorem 2.10 (Uniqueness). Let a a>t : [t,t + c) — >• X be a solution curve of an arc field $ 
with initial position a at time t, and let a btU : [u,u + c) X be a solution curve of $ with 
initial position b at time u. Then we have 

d{a at t{t + s),a b>u (u + s)) < e KAS d{a,b) + C\t - u\ a , for s G [0, c) 

where C is a constant depending only on c and 

Proof. Let us first define a a>T : [0, c) — >■ X by 

cr a , r (s) := o"q, )T (t + s), for s G [0, c) 

where a G X and r G [0, oo). If we can show 

d{?cM, °b,u( s )) < e KAS d(a, b) + C7|t - for s G [0, c) (2.27) 

for some constant C which is depending only on c and Ka then we are done. 
Triangle inequality gives 

d(a a , t (s), &b,u( s )) ^ d(a ajt {s), a bjt (s)) + d(a bjt {s), d- bjU (s)) (2.28) 

First, let us estimate the second term in the right hand side of (12.28j) . Let be the n-th 
discretized solution of a b j and let £^ be the n-th discretized solution of a bjU . We exploit 
Condition C and get 

di--=d(&(h,&{h)<£\t-u\° (2.29) 
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Again, triangle inequality gives 
4: = ^),^)) 

< ^„(|)^rr(^i)))+d($?r(^(i)),c(|)) 

= d(^¥r(^(^)),*?r(^(^))) +d(^r n (^)),$?r(^(^))) (2-30) 

We exploit Condition C to get 

d^rm^T-rm^))) < £\t-u\° (2.3i) 

and Condition A gives 

^rm±))) < d(^(i), ^(i)) (1 + |±) 

= *(l + fr) (2-32) 



We combine(E30]), ( l2~3H and ( l2~32l) 



z — 1 V on / On 1 1 



2n 1^ U \ \}~ 0- 2n ■ 



Similarly, for all k such that |^ < c, we have 

4: = 4ed),cd)) 



< - tt |«[l + (1 + ^1) + ... + + 
C , |t- M | Q [l + f^] fc - 1 



< 



2« 1 + *a _ 1 



C( e KA/c - 1) 

<^—- -±\t-u\ a (2.33) 

This means, for all s such that s < c 

d(e n (s),C(s)<C\t-u\ a (2.34) 
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where C := C(e KA ^ c — 1)/Ka- Since (12. 34[) is true for all n, we have 

d{o b>t (s),o b>u (s) < C\t - u\ a (2.35) 
Now, let us estimate the first term in the right hand side of (12.281) . We define 

g(s) := e- KAS d(a a , t (s),a bit (s)) 

For h > 0, we have 

g(s + h)- g(s) = e- KA ^d{a a , t {s + h), a b , t (s + h)) - e^d^s), o b , t (s)) 

< e - KA ^ +h \d(^ h +s (o a , t (s)), $^ +s (a M ( S )) + o(h)) - e- KAS d(a a , t (s), M*)) 

< e~ KAS e- KAh d(a a , t (a), cr M (a))(l + - e^ s d(^(s), cr M (a)) + o(/i) 
= (e" x ^(l + K A /i) - l)e- KAS d(a a , t (s), a b , t (s)) + o(/i) 

= o(h)(g(s) + 1) 

Hence, 

D^g(s) := hmsup( J < 

Consequently, g{t) < g(0) equivalently 

d(a a>t (s),a b>t (s)) <e KAS d(a,b), for s G [0, c) (2.36) 
(12T35D and (T2~36D together with d2~28l) give (12T271) and conclude the proof. □ 

2.3 Generating flows 

As a corollary of Theorem 12.91 and Theorem I2.10[ if an arc field $ satisfies Condition A,B 
and C then there is a unique solution curve a a ,t '■ [t,t + c att ) — > X with initial position a 
at time t, for each a £ X and t G [0, oo). To guarantee that c ajf = oo for all a G X and 
t G [0, oo), we borrow the idea of [I]. 

Definition 2.11. An arc field $ is said to have linear speed growth if there is a point x £ X 
positive constants ci(x) and C2(x) such that for all r > and t, I > 

p(x,t;r,l) < ct(x)r + c 2 (x) (2.37) 

Theorem 2.12. Let $ &e an arc field which has linear speed growth. Suppose that at each 
point a G X and t G [0, oo), $ has a solution curve o a ^ t '■ [t,t + c a> t) —> X with initial position 
a at time t. Then c aj4 can be chosen to be oo. 

Proof. Similar to Theorem 4.4 of jl] □ 



12 



3 Adding Flows 

3.1 Sum of arc fields 

Let $ and \1/ be two arc fields satisfying Condition A,B and C. We impose a certain com- 
mutation law on $ and 

Condition D: There exist constants < a < 1 and Cj > 0, and a function II : X x X x 
[0, 1] — > R such that for each a G X and £ G [0, oo), there are constants r a > 0, e a G (0, 1] 
and T t > t such that II is bounded from above on B(a, r a ) x B(a, r a ) x [0, e a ] and 

d(m h o n^&i), ° ^(62)) < w + hn(h, b 2 , h)) + c d h i+a (3.1) 

for all 61, 6 2 G -B(a, r a ), /i G [0, e a ] and s, s + h G [£, T t ]. 

Remark 3.1. Once we have fixed a G X, £ G [0, 00) and fixed constants r a ,e a ,T t then the 
function II is bounded above. We denote upper bound of II by K~d. 

Definition 3.2. For given a G X and £ G [0, 00), a solution curve of the sum of $ and \1/ 
with initial position a at time £ is a map a : [£, £ + c) — >■ X such that cr(£) = a and for each 
s G [£, £ + c) 

lim rfK g + 2/.),^+ fe o^(a(.))) _ Q 

fc->o 2h { ' 1 

Remark 3.3. If <3> and \l/ satisfy Condition D, then we can check that 

Hm d(<T(3 + 2h),9'+ h oV h Ms))) = lim dMa + 2h),& h +h onMs))) (3 3) 
/i->o 2/i h->o 2h 

So, a solution curve of the sum of $ and \1/ is a solution curve of the sum of \& and $. 

For a notational convenience in later computations, let us introduce new arc fields $ and 
\P which are defined by 

n(b) := n(b) := % h (b) 

for all b G X and (/i, s) G [0, 1] x [0, 00). It is trivial to see if <3> and \& satisfy Condition D 
then $ and ^ satisfy the following condition 

Condition D': For all 61, 62 £ r a) ; h G [0, e a ] and s, s + /i G [£, T t ], we have 

o §«(&!), o $*(& 2 )) < 6 2 )(1 + ^(61, 6 2 , /i)) + C^ 1+Q (3.4) 

with same constants r a , e a , T t , a, Cd and a function IT as in Condition D. 

It is also trivial to see that a : [£, £ + c) — >■ X is a solution curve of the sum of ^> and \l/ 
if and only if a satisfies 

[im «!(„(. + 2/,.), *r»o i-Ms))) 

h^o 2h v ; 

for all s G [£, £ + c). 
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Lemma 3.4. For a given a G X and t G [0,oo) ; let arc fields $ and ^ satisfy Condition 
A,B,C and D', and r a , e a and T t be constants in those conditions. If bi,b 2 G B(a,r a ), 
Ah G [0, e a ] and s, s + 4h G [t, T t ] then we have 

d(% +3h o o o *• fa), o ^(fc)) < dfa, 6 2 )(1 + />A) 3 + C(/i) 

w/iere A := max{A A , K B , K D ) and C(h) := C d h 1+a (l + AA) + hg(h,h)[l + (1 + AA) 2 ]. 




Figure 3: 

Proof. 1. We use Condition A to get, 

dt := d(<j> s h (h),<f> s h {b 2 )) < d(6i,62)(l + />A A ) 
2. Condition A and D give, 

< d(^ h (h), n(h))(l + hK D ) + C,/i 1+a 

< d(6i, 6a) (1 + hK A )(l + hKo) + C d /i 1+Q (3.6) 
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3. We combine Condition A and ( 13. 6p to get 

d 3 : = d« 3h (& h +2h o o $^( &1 )), n + 3h (n +2h ° n +h ° mm 

< d^ s h +2h o o S-fo), o o $^ 2 )))(1 + hK A ) 

< [d(b u b 2 )(l + hK A )(l + hK D ) + C d h 1+a ](l + hK A ) (3.7) 

4. Equation (12.6P gives, 

d A := d(^(6 2 ), $^ o < ^(fc, fc) (3.8) 

5. We exploit Lemma [2.41 to have 

d{n+3h q Q $ . (6a)))j ^(^(ft,))) < d ($j** o $*(& 2 )), $^(fo 2 ))(l + ^) 2 + ^(/i, h) 

and this together with ( 13. 8ft gives 

4 : = d(**+ 3ft o o v^($^ 2 ))) 

< hg(h,h)(l + hK A ) 2 + hg(h,h) (3.9) 

Finally, triangle inequality with ( 13. TP and ( 13. 9 ft gives 

° ^ ° ° ^(&l)> ° *3fc(k)) 

< d 3 + d 5 

< d(bi, + ^a) 2 (1 + hK D ) + C d /i 1+a (l + WT A ) + hg(h, h)[l + (1 + ^a) 2 ] 

□ 



Remark 3.5. By exactly the same argument as in Lemma 13.41 we can show a more general 
formula: Let h\ := h + r and h 2 := r + / 2 , then we have 

< d(6i, 6 3 )(1 + + ri^)(l + l 2 K A ) + C7(Zi, r, Z 2 ) (3.10) 

where C(h, r, / 2 ) := rg(h, r)(l + Z 2 i^)(l + rK A ) + Z 2 ^(r, / 2 ) + C d r 1+a (l + l 2 K A ). 

If we define K := maxji^^, K^} and u := max{/ 1; r, l 2 } then (13 . 1 0[) can be simplified as 

where C(u) := C d M 1+a (l + uif) + Ufif(u,u)[l + (1 + uK) 2 ]. 

3.2 Existence and Uniqueness of a solution for the sum of two arc 
fields 

Theorem 3.6. (Existence) Let <3>, \l/ be arc fields satisfying Condition A,B,C and D. Then, 
for given a G X and t G [0, oo), there is a solution curve a : [t,t + c) — >■ X of the sum of <3> 
and \& with initial position a at time t. 
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Proof. Without loss of generality, we assume $ and \& satisfy Condition A,B,C and D', and 
find a curve a satisfying (13. 5p . 

For a positive integer n, we define the n-th discretized solution by 



Us) 



* s-2 



*«(a) 

n (£n( 2" )) 



0< s< £ 



■i- < S < ± 

2 n — — 2" 



Qt+i-2 »/£ /J_NN J_ < g < 1+1 if 

^„ - n— n VS^V 2* 1 // 2" — — 2 n 



-i-2~ 



2m 



^ s t%Z(U^)) ^<s<^ if i = 2m+l 



Suppose r, I > are chosen so that p(a, t; r, /) < oo. If p(a, t; r, I) = 0, then a(s) := a defines 
a solution curve. Thus we assume p(a, t; r, I) > 0, and let 



c := mm 



p(a,t;r,l)' 



(3.12) 



It is easy to see that we have £ n (s) G -B(a, r) for < s < a This also implies {£, n }^ =1 is equi- 
Lipschitz with Lipschitz constant p(a,t;r,l). Moreover, by choosing r smaller if necessary, 
we may assume there are constants Ka, Kb and e G (0, 1] such that A(p, q, h) < and 
Q(p,l',h) < Kb for all p,q G B(a,r) and G [0, e]. We may also assume B(a,r) is a 
complete metric space. 

Let us first estimate the uniform distance between £ n and £ n _i. We apply Lemma [3.41 
with h := l/2 n , s:=t + §,b 1 := £ n (§) and b 2 := f n _i(|0 to have 

= o $^ o o A), nt 2h o n h (tn-i&)) 



Ai 
2~ r 



where K := max{K A , K D } and C(h) = C d h 1+a (l + hK) + fa/(7i, /i)[l + (1 + hKf] 
In general, we have 



< i(6.(%^),6^.(^5^)) (i + ^f 2 + c(l) [i + (i + 



k 2' 

/< 

2' 



< <K6.(o),6^i(o))(i + £) 3<i+1) + c(l) [i + (i + £) 3 + ■ • • + (i + £) : 
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Figure 4: 



Since £ n (0) = £„,_i(0), for all i such that 4(z + l)/2 n < c, we have 

^(^),e„- I (^))<c ( 2-»)( 1 + W- 1 



(1 + 2-"iT) 3 - 1 

< C(2~ n ) 



— K 1 

e * — 1 



2-«iT[(l + 2~ n K) 2 + (1 + 2"™K) + 1] 

< ^liH^l ,3,3) 

Notice, for sufficiently small h 

C{h) = C d h 1+a {l + hK) + hg(h, h)[l + (1 + hK) 2 } < 3h[C d h a + g(h, h)] (3.14) 
We assume n is large enough in ( I3.13p . We combine that with f )3.14p to show 

i(U^AUn-/^A)) < [C^r + K^,^)]A (3.15) 
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3c k 

where A := e ^ K ~ 1 ■ We notice A is independent of n. 
So for any s G [0, c), let i be an integer such that 

4z 4(z + l) 

— < s < — 

2" ~~ 2 n 

then we have 

4-i 4? 4? 4? 

<lp(^rJ)+[C4^T + H^^ (3-16) 

where we use the equi-Lipschitz property of £ n and (13.151) . 

Next we exploit (13.161) to show £ n is a Cauchy sequence in the uniform topology. For any 
s G [0, c), we have 

m— 1 



j=0 

8 1 a 1 1 

i=o 

8 °° 1 11 



j=n+l 

->■ as n -»• 



Since £ n (s) is in the complete space -B(a, r), we know that £ n converges uniformly. We define 
a : [0, c) — > X by the limit i.e 

a(s) = lim £„(s) 

n— too 

Next, we are going to show that 

to ■*(*(. + M), ^'°«r(*M)) 

for all s G [0,c). Once we have (13.171) . we are done with the proof by defining a solution 
curve a : [t, t + c) — > X as 

cr(s) := a(s — t), for sG [t,t + c) 

To prove (13.171) . we choose an arbitrary s G [0, c). Let e > and h > be fixed such that 
s + h < c. From triangle inequality, we have 
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d(a(s + 2h),*l +a+h o& h +8 (a(s))) < d(a(s + 2h),U^ + 2h)) , d(£ n (s + 2/i), o $* +s (£ n (s))) 



+ 



h h 



h 

Since £ n converges uniformly, we can choose n large enough so that 



(3.18) 



h + h "2 (3 ' 19) 

We combine ( 13~T8|) and ( 13~T9|) to get 

rf(a(s + 2/^) ! vI/^o^ +s (a( g ))) < rf(a( g + 2^),vl/^ +fe o$^(e ra (s))) + e 
/i ~~ /i 2 

We need to estimate the second term of (13.201) . Let I be a nonnegative integer such that 

21 ^ 2(1 + 1) 

< s < 



2 n ~ 2 r 



and define s' := J|, then 



d(Us + 2/i), tt t + s+ ' 1 o Q?' < d(£ n (* + 2/i), ^ + 2/i)) + d(^(a' + 2h),^+ s ' +h o $* +s '(^(0)) 

+ dp^-* o $r'(a( s ')), ^ +s ' +h o K +s '(Us))) 

+ d(^' +h o $^'(^(s)), 9?'+ h o 

< d(Us' + 2h), ^ S ' +H ° K +S '(Us'))) 
+ p(a,t;r,l)(s-s')[l + (l + hK A ) 2 } 
+ Ch(s- s') a [l + (l + hK A )} (3.21) 

where we use Condition C and the Lipschitz continuity of £ n in the second inequality. For n 
large enough, s — s' is sufficiently small so that 

p(a, t; r, l)(s - s')[l + (1 + hK A f] + Ch(s - + (1 + hK A )] e 

h " 2 { ] 

We combine ( 13T2TT) and ( T3~22l) to get 

d(en( g + 2/.),^ +fe o$^(^(s))) ^ w ( g / + 2^),vl/^ /+fe o«l>^ / (e w ( g 0)) 6 

" & + 2 ( j 
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Figure 5: 



We need to estimate the second term of (I3.23p . For the moment, let us assume h = 1/2 J 
for some j G N . With this assumption, we can apply Lemma 13.81 with m := n — j and 
b := £ n (s')j and § et 



n-j 



dW +h o K +S '(Us%Us' + 2/1)) < Ch[h° + p 



i=i 



C l a "~ J 1 



oo 



— 2,3 t ( 2^ ^ ~^ ^ ^ ^ ( 91 ' 91 ) ] 



i=3'+l 



2* 2* 



(3.24) 



Notice (I3.24p is independent of n, i.e it holds uniformly for large n E N. Now we combine 
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(J3HD, (EE) and and let ft(= l/2 j ) ->■ then 

[im jjgg + 2„), y ° »»(*(*))) < , im g[( 1 + f; ( 1 1 }] + t = t (3 25) 
h-+o ft ~j->oo LV 2-? ; Z-< y v 2 t'2t/J v ; 

This gives f[3~TTj) . 

For general ft, let be an integer satisfying 

— <h< — (3.26) 
2 n — 2 n 

and define ft' := We exploit Lemma [3.71 to estimate the last term in (I3.20p and get 
ft ~ h! 

- h> + 2 {6M) 

where we assumed n is large enough to get second inequality. We combine (I3.20p and f!3.27j) 

d(a(s + 2h),9 t +' +h o&+'(a{s))) ^ d(^ h t s+h ' o Us + 2ft')) | £ 

ft — ft' 

Let m be a nonnegative integer such that 2 m < k < 2' m+1 i.e 

2«i 2 m +i 

— < ft' < (3.29) 

On — 9 n 

From Remark 13.91 with b := £ n (s) and M := ft'//c, we have 



DC 



d(¥ h v +h, o^(Us)),us+2h'))<c2h'[(2hr+ E ( 3 - 3 °) 

i=m+2 

We combine (E2HD and (13T30]) to get 

^ + 2,),^o, r(gW )) £2c[(2fty+ £ i( i i )]+e (33i) 

i=m+2 

Notice that ft' converges to ft and m increase to oo as n — > oo. This implies, for all sufficiently 
large n, we have 

oo 

2C[(2ft') a + £ ^)] < 2C(2ft)« + 6 (3.32) 

i=m+2 
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which is independent of n. We combine (13.311) and (13.321) . and let h converge to then we 
have 



Um d( 9 (, + 2A),^c^(q(.))) < Um +e = e 

h^O h h->0 

This gives (I3.17P and concludes the proof. □ 

Lemma 3.7. Let £ n : [0, c) — > X be the n-th discretized solution constructed in Theorem \3.6i 
and < h < c be fixed. If k/2 n < h < (k + l)/2 n for some integer k then there is a constant 
C > such that 

d(n +h o n(Us)),Us + 2h)) d(Vt+ h 'o<S>UUs)),Us + 2h')) + c \h - h'\ 



h h h 



+ C\h-ti 

Here, h' = k/2 n and C depends only on c and Ka- 
Proof. Triangle inequality gives 

+ d{£ n {s + 2ti)^ n (s + 2h)) (3.33) 

Let us first estimate the last term of (13.331) 

d{£ n {s + 2h'),£ n {s + 2h)) < 2p(a,t;r,l)\h - h'\ (3.34) 

which comes from the Lipschitz continuity of £ n . 

To estimate the first term of (I3.33p . we exploit triangle inequality 

(3.35) 

Let us estimate the righthand side of (I3.35P term by term. For the first term, we use the 
Lipschitz continuity of $, \l/ and Condition A to get 

d(n +h ' o n(u*))> n^' ° mu*))) < d (n +h ' ° mu*))> n^ 1 ° mum 

< p(a, t; r, l)\h — h'\ 

+d^ u h {Us)),mus)m + hK A ) 

< p{a, t; r, l)\h - h'\[l + (1 + hK A )\ (3.36) 



22 



For the second term, we use Condition C and get 

din+h $£(£„(*)), n+ v mn{s))) < Ch{h _ hT (3 37) 

Combine fl3T35D .f l336ll and (EOT!) 

d(^i +h o n(U*)), n? h ' ° < *; r, 01* - h'\ [i + (i + 

+ (3.38) 

Finally, fl3T33|) . f lBTMj) and fl3T39|) give 

d(tf«+* o $«(£„(*)), £„(* + 2h)) < d(^+ h ' o (Us)), Us + 2/1 1 )) 

+ p(a, t; r, Z)|/i - /i'|[3 + (1 + Wi A )] + C7i|/i - h'\ a (3.39) 

Which concludes the proof. 

□ 

Lemma 3.8. Let $ and \1/ 6e arc fields satisfying Condition A,B,C and D'. For given a G X 
andt G [0, oo), letr a ,e a andT t be constants in those conditions. Forb G B(a,r a ), h G [0, e ] 
and s,s + h G [£, It], i/iere exists a constant C > suc/i £/&at 

1=1 

where m and n are nonnegative integers satisfying m < n. Here, C depends only on K and 
T t . 

Proof. We use Lemma [3. 4[ 

dl := d(^ h +h o n(b), ° *V? ° KT ° *W 6 » < C(\) 

where C{h) := C d h 1+a {l + hK) + hg(h, h)[l + (1 + hK) 2 ]. 

d 2 : = d^ h f 2 ° K% h o KT ° *W&)« <* 22 ° ° • • • ° ^ ° 

<^(l + ^ 3 + C-(A) (3.40) 

where J 2 := d^T ° ^/#), ° ° ^ ° ^(6)). 

Again, by Lemma 13.41 

d2 < C(A) (3.41) 
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Figure 6: 



We combine QS3DD and (jOIjl to get 



^ 2 <^)[1 + (1 + ^) 3 ] (3-42) 

Similarly, 

^ [i + (i + |*) 3 + (i + + (i + |*f 3 ] 
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In general, 



< C(£) [! + (!+ + (1 + + ■•• + (! + ^A')' 2 '- 113 ] 



Notice C(h) < 3[Cdh 1+a + hg(h, h)} for small h. And there is a constant C such that, for all 
I 



So we have 



This implies 



*<G![( 5 ) + 2f9(2f. y)] 



E*2cti[(^) 1+ °+|44)] 

i=l i=l 

i=l t=l 

8=1 



□ 



Remark 3.9. Lemma [3.81 says, if fc = 2 m then 



° ^ ab - 1 )« o . . . o ^ o < c{ku) [{kuY + £ sfe §)] 



2 4 2 l 

i=l 



In general, we can show that if 2 m < fc < 2 m+1 then 

<W" ° K H2k ' 1)U o o.-o^o $*(&)) 

< C(2^) [(2™* 1 ,*)- + £ $(1— ^, ±-J±)] (3.43) 



i=i 
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Theorem 3.10. (Uniqueness) Let a a)t '■ [t,t + c) — >■ X be a solution curve of the sum of two 
arc fields $ and \& with initial position a at time t, and let a^ u : [u,u + c) — >■ X be a solution 
curve with initial position b at time u. Then we have 

d(a a , t (t + s),a biU (u + s)) < e KAS d(a,b) + C\t-u\ a , for s G [0, c) 

where C is a constant depending only on c and Ka- 

Proof. Similar to the proof of Theorem 12.101 □ 
3.3 Adding flows 

Like what we did for a single arc field in section 2, we impose the linear speed growth 
condition on $ and and get the following theorem. 

Theorem 3.11. Let $ and \& be arc fields with linear speed growth. Suppose that at each 
point a G X and t G [0, oo), there is a solution curve of sum of $ and ^ 

a ajt : [t,t + c a)t ) -> X 

with initial position a at time t. Then c a ^ t can be chosen to be oo. 

Let A be the set of all time dependent arc fields which satisfy Condition A,B,C and linear 
speed growth condition. For each $ G A, solution curves of $ generate a time dependent 
flow and we denote it by er$. Likewise, we use notation cr$ + # for the flow generated by the 
solution curves of the sum of $ and when they satisfy Condition D. Notice that we have 
a^ + q, = er*+$ by symmetry. 

Let us define an equivalence relation ~ in A as follows 

$ ~ $ if cr$ = o~$ 

and we denote the equivalence class containing $ by [<E>], i.e 

[$]:={$ G A : $ ~ $} 

It is easy to see that cr$ can serve as an arc field and cr$ G [$]. 

From the argument above, there is an one to one correspondence between A/ ~ and the 
set of flows satisfying Condition A,B,C and the linear growth condition. 

T := {o\j> : $ G A} — {[$] : $ G A} = A/ ~ 

Definition 3.12. Let $, ^ G .4. and suppose that $ and \1/ satisfy Condition D. We define 
0$ + <T^ := 0"$+* and call it the sum of two flows cr$ and 0\p. 

Lemma 3.13. T/ie sum o/ £u>o /?ou>s i.e er$ + er^ we// defined. 
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Proof. Let $ 6 [<&], ^ 6 and suppose that $ and \I/ satisfy Condition D. We need to 
check 

0$+* = cr$ + ^ (3.44) 
To show (I3.44p . it is enough to prove 

lim d(n +h on(b),n +h on(b)) = o 

/i->0 /i 

for all 6 G X and s > 0. 

Let o"^, : [s, oo) — > X be the solution curve of arc field \1/ with initial position b at time s. 
From Condition A, we have 

Um d(^ h +h (n(b)),K +h (Ms + h))) < ^ d(^(6), + h))(l + hK A ) 



h^>-0 h h 

= (3.46) 



Similarly, we have 



h^O h /wO h 

= (3.47) 
We combine (ETIB]) . fl3T47j) and = <7$, [$] = [$] to get 

ft.— s>o h ~ h^o h 

= 

which gives (I3.45P and concludes proof. □ 

Now, let us think about three flows cr$,a$ and ere. Suppose $ and \l/ satisfy Condition 
D, then the flow 0"$ + cr#(= cr$ + ^) is well defined. Furthermore, if $ + ^ and satisfy 
Condition D, then 

(T($ + ^) + e := (cr$ + er^) + <Te = er$+# + 0"e 
is also well defined. It is not hard to see that 

lim ggg+gneCf + 3fe), e-r ° ° n h ns))) = f3 . 

ft->o 3/i 
Similarly, if \1/ and satisfy Condition D, and $ and ^ + satisfy Condition D then 

o"$+(*+e) := 0$ + (c* + o-q) = cr$ + cr^ +e 
is well defined. We also have 

Hm g^Wj + 3fe), 0-+ 2fe ° m h ° gfcKf))) = f 3 491 

h^o 3h 1 ' ' 

By combining (I3.48P and (I3.49p . we have the following associative law in the sum of flows. 
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Corollary 3.14. Let $,$,0 6 A, if they satisfy properties related to Condition D to make 
sense of sum of three flows, then we have 

cr$ + (o> + er e ) = (cr$ + cr^) + cr e 
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